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By comparing the voltage-current (V — I) curves obtained before and after cutting a sample of 2H- 
NbSe2, we separate the bulk and edge contributions to the transport current at various dissipation 
levels and derive their respective V — I curves and critical currents. We find that the edge contribu- 
tion is thermally activated across a current dependent surface barrier. By contrast the bulk V — I 
curves are linear, as expected from the free flux flow model. The relative importance of bulk and 
edge contributions is found to depend on dissipation level and sample dimensions. We further show 
that the peak effect is a sharp bulk phenomenon and that it is broadened by the edge contribution. 

PACS numbers: 74.25.Fy, 74.60.Ge, 74.60.Jg, 74.60.Ec 



I. INTRODUCTION 

One of the remarkable differences between a normal 
metal and a superconductor stems from the way they 
carry current. While in metallic samples the current den- 
sity is usually homogeneous throughout the material, in 
superconductors in the Meissner state the current flows 
along the sample surface and edges so as to eliminate the 
self-induced magnetic field within the sample volume. In 
type-II superconductors in the mixed state, where the 
material is permeated by magnetic flux lines (vortices), 
a similar current enhancement along sample edges is the 
result of a surface barrier which inhibits the entry or 
exit of vortices p[-fi3|. Several types of surface barri- 
ers have been identified including the Bean-Livingston 
jjj and geometrical barriers. The former, which is 
the primary source of edge currents in the experiments 
described here, is a result of the competition between 
the attraction of a vortex to its image and the repulsion 
arising from its interaction with shielding currents. Re- 
cent Hall probe measurements in Bi 2 Sr 2 CaCu208 @-§] 
and NbSe 2 crystals have shown edge current enhance- 
ment due to the surface barrier. But thus far these edge 
currents were only qualitatively identified. 

The experiments described here allow for the first time 
to derive individual voltage-current (V — I) character- 
istics of the edge and bulk currents. Below the peak 
effect region (a peak in the critical current just below 
T c ) we find that vortex entry and exit at the edges is 
governed by thermal activation across a current depen- 
dent surface barrier, whereas vortex motion in the bulk 
is non- activated. Our results show that the observed 
nonlinearities of the V — I characteristics are due to the 
edge contribution while the bulk V — I curves are lin- 
ear. We further show that the peak in critical current 
is primarily a bulk effect which sharpens and becomes 
more pronounced when the edge contribution is removed. 



These experiments demonstrate that boundaries have a 
profound effect on the V — I characteristics as well as 
on the field and temperature dependence of the critical 
current. 

II. EXPERIMENTAL DETAILS 

The experiments were carried out in the low tem- 
perature superconductor 2H-NbSe2 where local self-field 
and magnetization |Dj ] measurements unveiled the 
presence of significant edge currents. Transport, neu- 
tron scattering and magnetization measurements in these 
samples [[l4]-^6| revealed a number of unusual phenom- 
ena including memory and current driven reorganization 
]l(i| plj . In addition, the shape of the peak effect in this 
material was found to change significantly with measure- 
ment speed or contact configuration [ p5[ . Sev- 
eral of these results were shown to be a consequence of 
vortices traversing a surface barrier as they enter the 
sample |2jj2^]. In the experiments described here the 
surface barrier was determined after separating the edge 
and bulk contributions to the current. The separation 
procedure involves cutting a sample to reduce its width 
and comparing the V — I characteristics before and after 
cutting. A schematic illustration of the cutting is given 
in Fig. 1(a). 

The data were acquired on two undoped single crys- 
tals of 2H-NbSe 2 with initial sizes of 8(L) x 1.72(w) 
x 0.020(d) mm 3 (sample A) and 6.3(£) x 1.40(w) x 
0.060(d) mm 3 (sample B) and with zero field critical 
temperatures and width of T c = 7.18 K, AT C = 95 mK 
and 7.21 K, 92 mK, respectively. Our measurements em- 
ployed a standard four probe technique with low resis- 
tance Ago.iIn .9 solder contacts. The distance between 
the voltage contacts was 2.5 mm and 1.5 mm in samples 
A and B, respectively. The critical current, I c , is defined 
as the current at which the voltage reaches 1 //V. The 
magnetic field was kept along the c axis of the sample and 
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FIG. 1. Schematic illustrations of sample cutting (a) and 
the corresponding current distributions in the vortex state 
(b).The shadowed area represents the edge current I 3 (see 
text). The inset in (b) shows the resistance versus tempera- 
ture curves and the determination of the ratio of the sample 
width before (wo) and after (wi) cutting. 



the dc current was applied in the a — b plane. The vor- 
tex lattices are prepared by applying the magnetic field 
after cooling the sample through T c (zero- field cooling). 
Sample cut was carried out manually with a sharp razor 
blade. The width reduction factor ao n = wo/w n (wo, 
w n are the width before and after n th cutting) was de- 
termined by direct inspection under a microscope and 
confirmed by measuring the ratio of normal state resis- 
tances before and after cutting. Sample A was cut twice 
with aoi = 1-52 [shown in the inset of Fig. 1(b)] and 
<^02 = 3.65 for the first and second cuts respectively, and 
sample B was cut once with «oi — 1-35. By using this 
procedure rather than samples with different widths one 
can be certain that in comparing the V — I curves before 
and after cutting all the parameters (excepting the newly 
cut edge) are the same. 

III. RESULTS AND DISCUSSION 
In Fig. 2 we show the effect of sample cutting by com- 
paring the V — I curves and critical currents before and 
after cutting. In Fig. 2(a) the results are shown for 
both the normal (T = 7.6 K > T c ) and superconduct- 
ing (T = 4.25 K < T c ) states at H = 1 T. A compar- 
ison of the measurements at the same average current 
density (J = I/dw) is obtained by plotting the voltage 
against the scaled current «02^ (solid line). At 7.6 K 
the scaled curve exactly overlaps the response in the un- 
cut sample, clearly showing that the current density in 
the normal state scales with the inverse of the sample 
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FIG. 2. V — I curves at temperatures below (4.25 K) and 
above (7.60 K) the critical temperature T c (a) and temper- 
ature dependence of critical current (b) for sample A before 
and after sample cutting. The open symbols and solid lines 
represent the as-measured and the scaled data. 



width. In other words the current distribution in the 
normal state is uniform. Using the same procedure for 
the data in the vortex state, we find that the scaled V — I 
curve of the cut sample is shifted to the right of the ini- 
tial curve, indicating that for the same voltage response 
the average current density in the cut sample is much 
higher. The same tendency is found for the scaled criti- 
cal currents, which are significantly larger after cutting, 
Fig. 2(b). As shown below these results are consistent 
with an enhanced current at the sample edges due to 
the surface barrier. In fact both theories [p] JlO|Jll[ | and 
experiments [^-[7 p2|Jl~^ ] favor an excess current carrying 
capacity at the sample edges. The results are analyzed 
by separating the total current, I = I s + lb, into a uni- 
formly distributed bulk contribution, lb = Jbwd, with 
current density Jb and a nonuniform contribution, I s . 
The latter, assumed to remain unchanged after cutting, 
represents the edge current due to the surface barrier 
|l|-|| . If edge contamination |24|] is present the resulting 
enhanced edge current would also be included in I s . This 
model is applicable when the sample width after cutting 
is much larger than the characteristic extent of edge cur- 
rents. In our experiments this condition is well satisfied 
below the peak effect regime where edge currents are pre- 
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FIG. 3. V — I curves for total (thick line), edge (dashed) 
and bulk (thin line) current before (a) and after (b) cutting 
the sample, (c) edge V — I s curves in 1 T and at various 
temperatures 



dominantly due to the surface barrier (edge contamina- 
tion is negligible p4|j ) an d are thus confined to within 
a narrow strip < A ~ lOOnm <C w [[10 11|. It follows 
that the vortex velocity and hence the voltage response 
at a given field and temperature is uniquely determined 
by the bulk current density Jf, and by I s . In Fig. 1(b) 
we plot a schematic current distribution in the presence 
of a surface barrier [jl 10 ll[] or/and edge contamination 
illustrating the effect of cutting for a constant Jf,. In this 
model the voltage response to a driving current / = It+I s 
in the initial sample will be the same as the response to 
a current I a = It/ a + I s in the cut sample, leading to 
a straightforward procedure for analyzing the data and 
for separating bulk from surface currents. Thus, after 
obtaining the values of / and I a at a given voltage re- 
sponse by measuring the V — I curves before and after 
cutting the sample, the bulk and surface contributions 
to the current are given by: h — a(I — I a )/(a — 1) and 
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FIG. 4. Temperature dependence of total, edge and bulk 
critical currents at 1 T. Main panel: results for sample A 
(u>2 = 0.471 mm); inset: expanded view of the peak region 
for sample B (wo = 1.40 mm). 



I s = (al a —I)/(a — l). Repeating the same procedure for 
various voltage levels, separate V — I characteristics can 
be obtained for the bulk (V — lb) and the edge (V — I s ). 
Results for the V — I characteristics of sample A at 4.25 
K and 1 T are shown in Fig. 3 (a) and (b) for the initial 
(wq = 1.72 mm) and cut sample (w2 = 0.471 mm) respec- 
tively. We note that the relative contribution of the edge 
is larger in the narrower sample Fig. 3(b) than in the 
wider one Fig. 3(a), in accord with the assumption that 
I s is unchanged, whereas lb is proportional to the sample 
width. In addition, the relative contribution of the edge 
to the total transport current diminishes with increasing 
dissipation level indicating that the voltage grows much 
faster with increasing I s than it does with lb- We fur- 
ther note that the measured V — I curves in both the 
initial and cut samples are nonlinear at modest dissipa- 
tion levels. Previous interpretations of the nonlinearities 
in the V — I curves in 2H-NbSe2 were usually based on 
the assumption of a bulk phenomenon. However from 
Fig. 3(a) and (b), our data show that the nonlinearity 
is due to the edge contribution whereas the bulk V — I 
curves are remarkably linear, consistent with early mod- 
els of vortex motion in low temperature superconductors 



Theoretical studies of vortex entry and exit across sam- 
ple edges |l(],[ll| have shown that a surface barrier U(I S ), 
gives rise to thermally activated V — I s characteristics: 
V = Voexp[—U(I s )/kBT]. The current dependence of 
the surface barrier is determined by the mechanism of 
vortex penetration. For example, the surface barrier in 
a slab geometry was found to have a power law current 
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dependence, U ~ 7 S , in the 3D case while in the 2D 
case it is logarithmic U ~ U ln(I / Is) 1 / 2 with f/ a char- 
acteristic energy scale In Fig. 3(c) we plot the 
edge V — I s characteristics at 1 T for several tempera- 
tures. The data are best fitted with a logarithmic current 
dependence of the surface barrier U(I S ), as seen from the 
straight lines obtained on a log-log scale. The tempera- 
ture dependence of the slope is consistent with thermal 
activation with U ~ 70 K for T < 5 K and decreasing 
in value for T > 5 K. 

The temperature dependence of the critical current for 
samples A and B is shown in Fig. 4. The as-measured 
critical currents are shown together with the bulk (I c b) 
and edge (J cs ) contributions obtained by the procedure 
described above. The two samples were grown in the 
same batch and are expected to be similar in quality, 
but they differ in thickness, with sample B three times 
thicker than sample A. Comparing the bulk critical cur- 
rent densities (J cb = I c b/wd) in the two samples we find 
that below the peak effect they are practically identical, 
despite the fact that the average critical current density, 
Jc = I c /wd, calculated in the usual way as the total crit- 
ical current divided by the sample cross section is more 
than 3 times larger in sample B than in A. These data 
show that the two samples are identical in their bulk 
pinning properties and that the distribution of pinning 
centers in the bulk is homogeneous. The significant dif- 
ference between the bulk and average critical current den- 
sities must therefore be due to the surface barrier. Com- 
paring the surface critical sheet current densities, defined 
as J cs = I cs /d, we find a strong enhancement of the edge 
critical current density in sample B compared to sample 
A, which implies that the surface barrier in the thicker 
sample is larger. More work is needed to understand the 
dependence of the surface barrier on sample thickness. 

We next consider the peak effect region [Q, where 
both samples exhibit a well defined enhancement in the 
total critical current by factors of 2.0 and 2.7 for sample 
A and B respectively. The shape of the peak was previ- 
ously found to depend on measurement speed and contact 
geometry becoming much sharper and more pronounced 
when the contribution from vortices entering through the 
edges was reduced or eliminated |24|,|25| . A similar result 
is obtained here by separating the edge and bulk contri- 
butions to the critical current. As shown in the inset of 
Fig. 4 the bulk peak effect consists of a 20 fold enhance- 
ment in I cb [I cb {T = 5.63 K)= 2.9 mA; I cb (T = 5.66 
K)= 59.6 mA] over a temperature range narrower than 
the width of the zero field superconducting transition. 
This remarkably sharp peak is similar to that obtained 
when vortex crossing of edges is eliminated by using a 
Corbino geometry |25|]. In the lower part of the peak 
I cs increases rapidly with temperature compared to its 
monotonic decrease below the peak. This enhancement is 
consistent with the edge-contamination mechanism sug- 



gested by Paltiel et al. ]23] since the higher critical cur- 
rent of the disordered phase, injected at the sample edge 
in the lower part of the peak effect, leads to an addi- 
tional contribution to I cs over and above that due to the 
surface barrier. These results are also consistent with 
the significant sharpening of the peak effect observed in 
high frequency ac measurements Jl8| , p~9|j2^ ] where edge 
contamination is practically eliminated. 

In the above analysis we assumed that the edge contri- 
bution is unchanged by cutting the sample. As a check 
we cut the sample a second time and compared the bulk 
V — I characteristics to those before cutting. Below the 
peak region we found that the V — J b curves were un- 
changed within better than 10%. 

IV. CONCLUSIONS 

In conclusion, by comparing the V — I curves of sam- 
ples before and after they are cut, we obtained the bulk 
and edge contributions to the transport current at vari- 
ous dissipation levels. This led to the first derivation of 
separate V — I curves for the bulk and the edge. Below 
the peak region the nonlinearity of the measured V — I 
curves in 2H-NbSe2 is due to the edge contribution and 
the edge current is governed by thermally activated vor- 
tex crossing through a current dependent surface barrier. 
By contrast the bulk V — I characteristics are linear con- 
firming the free flux flow model for vortex motion in the 
bulk. In the peak effect region the temperature depen- 
dence of the bulk critical current exhibits a very sharp 
peak. The edge contribution starts increasing before the 
bulk current does, leading to a smeared out peak in the 
total critical current. 
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